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Abstract 

Complex flavour couplings (off-diagonal mass terms) in the squark sector of super- 
symmetric theories may drastically alter both the rate and the CP-violating asymmetry 
of certain S-meson decays. We consider the effects of couplings that induce b — > s tran- 
sitions and lead to final state with strangeness one. We investigate the bounds that 
must be satisfied by the new terms and explore the possible implications on direct and 
mixing induced asymmetries in the charged and neutral B — > J/^K and B — > <pK 
decays. 
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1 Introduction 



Certain decays of P-mesons are expected to shed light on the mechanism of CP violation and, 
more generally, on new physics. Indeed, the first surprising results on the CP-asymmetry 
a^x in the decay B — ► J/^/K raised the hope for a first new physics signal. Meanwhile, the 
value of the measured asymmetry has changed considerably and the present world average IT] 

a J™ = 0.79 ±0.12 (1) 

almost coincides with the standard model expectation |2| 

~ 0.70 ±0.10 . (2) 

Although this came as a disappointment, it reinforced the view that one should be prepared, 
both on the theoretical and the experimental side, for unexpected observations. In fact, the 
small error of the above world average and the improvements that are expected during the 
next years will allow to detect even small deviations from the standard model predictions. 

There are of course many new physics scenarios. On the one hand, it is possible to 
describe their signatures in very general terms (for recent expositions, see e.g. Refs. P-H). 
On the other hand, one can propose a specific new physics model and investigate its con- 
sequences. The latter approach was the topic of a huge body of work. For what concerns 
supersymmetric models, emphasis was given recently to the so-called minimal flavour vio- 
lating models (MFV) and their possible variants |jj . 

In virtually all new physics models, new non-standard fields are introduced and, there- 
fore, new complex (i.e. with CP-violating phases) couplings appear. This is well known 
in supersymmetry to which we turn for definiteness. In supersymmetric models there are 
several classes of phases. Those in the \x and flavour diagonal A terms do not contribute 
to flavour changing processes but are strongly bound by the electric dipole moment of the 
neutron and other particles ||. The phases of the Yukawa couplings are the same as in the 
standard model. If these are the only phases and flavour changing couplings of a model, 
the latter belongs to the class of MFV supersymmetric models that exhibit many analogies 
to the standard model ||. Therefore, the most interesting phases are those in the squark 
and slepton mass matrices. Their flavour diagonal elements are either real by definition (in 
the LL and RR sectors) or small (in the LR and RL ones) as pointed out before (restric- 
tions from electric dipole moments). Therefore the flavour changing elements (i.e. {m 2 )^ 
and (m 2 )if, in the sectors LL, RL, etc.) are the most interesting ones and we will focus 
on their effects. They contribute mainly through loop diagrams with internal gluinos and 
charginos. Since we are interested in contributions to Wilson coefficients that are already 
quite large in the SM |(i.e. the coefficients of the QCD penguin and chromo-magnetic dipole 
moment operators), we will only discuss gluino loops; in fact, they tend to dominate over 
the corresponding chargino ones. 

* Contributions to the small electroweak penguin coefficients are usually negligible because they do not 
appreciably change the standard model predictions. 
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Clearly this is not the first analysis of such SUSY diagrams. Their impact on flavour 
physics is known since a long time JJ, and there are numerous recent investigations. To 
clarify the new aspects provided in this paper, we review briefly some of the most recent 
efforts. 



In Ref. [10] a first comprehensive analysis of the gluino exchanges was given. The sys- 



tematics of the perturbative expansion was investigated in Ref. [11|, which also included an 



analysis of the mass insertion approximation: the latter was found to be sufficient in many 



cases and we will use it again in this study. In Ref. |L2[ the bounds from the rare decay 
b — > s-y were derived; however, the various couplings were assumed to be real. In Ref. |13[], the 
influence on the electromagnetic penguins (isospin violating terms) and the consequences on 
B — ► Ktt decays were studied; particular focus was given to the determination of the CKM 
angle 7. In Ref. || effects of new terms on direct CP- violating asymmetries in charged decays 
were given, but the possible values of the new coefficients were not investigated. In Ref. 
a completely general parameterization of new physics effects in the decays B — > <pK and 



B — > J /^K was given; however, no particular model was explicitely studied. In Ref. [14 
the case of left-right symmetric models were investigated and it was found that, in these 
models, sizeable deviations from the relation a^ K = cl^k are possible. 

In this paper, we consider the complex (2, 3) entries in the down squark mass matrix 
and work out their consequences for the time dependent CP asymmetries in the decays 
B — > <pK and B — ► J/^K. In the standard model, these asymmetries are equal and are a 
measure of the phase of the B^-B^ mixing amplitude (which in the SM is 2/3). The complex 
couplings that we consider leave the Bd-Bd mixing phase unchanged but have a significant 
impact on the amplitudes of the decays. The B — > J/^/K amplitude is dominated by SM 
tree level contributions and, consequently, the CP-violating asymmetry is hardly modified; 
on the other hand, there is room for such an effect in the transition B — ► <pK. This decay 
is easily accessible at the P-factory experiments [|l5j and an investigation of the related CP 
asymmetries is clearly worthwhile. In the analysis we will concentrate on this decay. 

This note is organized as follows. In Sec. |2]we review the effective Hamiltonian for b — > sss 
transitions in generic extensions of the standard model and give the explicit expressions of 
the various CP asymmetries. We also discuss the observables BiB — > Xg'y) and AMb s which 
are related to the b — > s transition and provide interesting additional pieces of information. 
In Sec. we present the SUSY model that we consider and the explicit contributions to the 
Wilson coefficients. The numerical analysis is presented in Sec. f|. The impact of the complex 
SUSY parameters on the branching ratio and CP asymmetry of B —>■ X s ^ is explored and 
we analyze the correlation between the mixing-induced (and direct) CP asymmetries in 
B — > (j)K and the B s -B s mass difference. A brief summary of our results and some comments 
are found in Sec. |5|. 
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2 Effective Hamiltonian 



The effective Hamiltonian for b — > s transitions can be written as 

10 



71 
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where (e U)C = 2/3 and e^. 
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In Eq. (0), we did not write the operators C" 2 obtained by the replacements c — > u in 01,2, 
and the semileptonic current-current operators which induce the transitions b — > s£ + £~ . The 
effective Hamiltonian for b — > d transitions can be obtained via the substitution s — > d. Also, 
operators with different helicity structures are not explicitely written and can be obtained 
from the above operator basis via the replacement L <->• R ||11|| . We will comment on their 
impact on the numerical analysis in section f|. 

The SM values of the Wilson coefficients Ci_io(m&) and C-j 1 ^gij n b) are given in table [I]; 
the opposite chirality operators do not get standard model contributions if the light masses 
are neglected. 



Ci 


-0.171 


CV 


-0.00001 


c 2 


1.070 




0.0005 


c 3 


0.0114 


c 9 


-0.01 


c 4 


-0.0321 


Cio 


0.0002 


c 5 


0.00925 




-0.313 


c 6 


-0.0383 




-0.188 



Table 1: SM Wilson coefficients at the scale nib. 



We will use the following approximate form | for the B — > (J/*&, (f))K amplitudes ||: 

2 

A* = -V2G F UF^ K {m%)m i) {e^-p B )V tb V* + ^^f C 8g r 8 ^ , (16) 



47r g 



Here, is the decay constant of the J/\l/(0), Ff^ K (q 2 ) is the B ^ K penguin form 

factor, Sumk — —0.76 is the ratio of the chromo-magnetic to penguin form factors for the 
B — > (J/^f,cj))K decay, q 2 ~ ml/2, and r 8 ~ 1/12 is the ratio of colour octet and singlet 
matrix elements. P is an 0(a s ) contribution to the matrix elements of the QCD penguin 
operators, it is insensitive to new short distance physics, and it is originated by loop diagrams 
with an internal charm quark. Since the typical q 2 is above the charm production threshold, 
P carries a strong phase and its numerical value is —0.0132 — i 0.0145 [|ll|. Finally, the 
coefficients are 



= Ci + C 3 + C 5 + ° 2 + ^ 4 + Cfi + 2r 8 (C 2 + C 4 + C 6 ) 
n r j r L p i ^ + ^4 + ^ 1 / „ „ „ C 8 + C 9 + C 10 

<^4> — <-/ 3 + 4 + 5 H I O7 + O9 + O10 H 



(18) 
(19) 



where all the C, are to be evaluated at the scale ra&. 

New physics contribute to the above Wilson coefficients; the resulting new phases will 
show up in the direct and mixing induced CP asymmetries. Before writing an explicit model 
for such contributions, we briefly discuss the renormalization group (RG) running of C^^fj) 
from a = 0(Mw) to /i = 0(rrib). The solution of the NLO RG equations reads: 

C^(m b ) = 0.3515 (1 + 0.0018 i? 3 - 0.0014 i? 4 + 0.0020 i? 5 + 0.0053 i? 6 ) , (20) 
C^nib) = -0.0243 (1 - 0.030 R 3 + 0.085 R 4 - 0.028 R 5 + 0.044 R 6 ) (21) 

where Ri = Ci(M w ) / Cf M \M W ) and the other coefficients are fixed to their SM value. It 
follows that the impact of the QCD penguin matching conditions (the values of the Wilson 

§This expression is based on naive factorization and is therefore not exact; but it correctly describes the 
way new physics enters. 
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coefficients at the scale M w ) on C^{m b ) and C^uib) is respectively of order 1% and 10%. 
Thus, as it is well known, their effects on C^{m^) is totally negligible: the coefficients of the 
current-current operators totally dominate the amplitude even if their contribution is colour 
suppressed. Thus, the SM prediction of an extremely small direct CP asymmetry in the 
J mode is unaltered by new physics contributions. On the other hand, new physics whose 
contributions are large enough to dominate the phases of C^-q{M\v)i may modify the overall 
phase of C^m^) at the 10% level. Taking into account that the CP asymmetry depends 
on twice the phase of C^rrib), we see that in principle large (0(0.2)) deviations from the 
SM relation a^ K = a^x are possible. Moreover, the presence of a term proportional to C 8g 
in Eq. (|P7j) can have a strong impact if there are large new complex contributions to the 
chromo-magnetic dipole operator. In the standard model, this term is negligible. 

Let us introduce the time-dependent CP asymmetry in the decays B d /B d —>■ <f)K, given 

by 

F(B° d (t) -> <j>K) + T(B° d (t) ^<f)K) CP V d ' CP V d ' V ' 

where AM^ is the is the Bd -Bd mass difference and B d (t) (B d (t)) is the state at time t 
which started as a pure B® {B®) at t = 0. Aq P and Aqp are the direct and mixing-induced 
CP-asymmetries, respectively. Their explicit expressions are 



-4cp = (23) 



21mA 



1+ A 



<f>K 



with 



(25) 



A 

A ' 

Here, (3 is the inner angle of the unitarity triangle of the standard model, 9d is a possible new 
physics contribution to the phase of the Bd — Bd oscillations, and A (A) are the amplitudes 
of the decay in question (and its CP conjugate). 

We turn next to the branching ratio and CP asymmetry in the decay B — > X s 7 and to 
the B s -B s mass difference. These observables are strongly affected by new physics in the 
FCNC b — >• s transition and therefore must be properly accounted for. 



2.1 B -> X s7 

The inclusive transition B — > X s 7 plays a major role in limiting possible new physics con- 
tributions to B decays both through constraints on the branching ratio and on direct CP 
asymmetry. The experimental information on the latter quantities is [|T7| 

B(B -> X sl ) = (3.22 ± 0.40) x 10" 4 (26) 
A CP {B^X sl ) = (-3.5 ±7.7)%. (27) 
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In the numerical analysis we will use the NLO computation of B(B — > X^) which, as a 
function of the Wilson coefficients evaluated at the scale M\y reads 



B = [1.258 + 0.382\R 7 \ 2 + 0.015\R 8 \ 2 + 1.395 ReR 7 + 0.161 ReR 8 + 0.083 Re(R 7 R* 8 )] 10~ 4 

(28) 

where R 7t $ = C 7 ^ Sg (Mw) / C 7 ^f 8g (Mw) ■ The parameter 6 that determines the cut on the 
photon energy spectrum is set to 0.90 according to Ref. |18| . 

Similarly, the CP can be written as fi9|| : 

A CP = 1.06 Im— - — r - 9.52 Im— *j — - + 0.16 Im — — \ . (29) 



2.2 AM B 



The effective Hamiltonian for the AS = 2 transitions can be written as 

2 



n/AS=2 
H eff 



where 



GlMl 



1 w iv tb v*T 



qVllqvll _j_ ^qSllqSll _|_ qSlrqSlr^ 



i=i 



+ h.c. 



(30) 
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Q SLL 
qSLR 
SLR 



[s L ^b L )(s L -f^b L ) 
{s L b R )(s L b R ) , 



[s L b R )(s R bi 



R) 



(31) 
(32) 
(33) 
(34) 
(35) 
(36) 



together with the operators O yRR and Of RR , obtained from the corresponding LL ones via 
the substitution L <-> R. 

In order to minimize the impact of hadronic uncertainties we will consider, as usual, the 



ratio X s d = AM^/AMg^. Its explicit expression is [|20] 



X 



sd 



Mb, 



V, 



ts 



Vt 



td 



c s 
c 1 



where 



C s = C 



VLL 



OKf. 



qSLL _j_ ^±(jSLL 



K K 



c SLR + ^ Ks - c SLR 



(37) 



(38) 



with n s = Mg /(m b + m s ) 2 and £ is a ratio of hadronic matrix elements, numerically equal 
to 1.16 ± 0.05 ||. The coefficient C d is obtained by replacing s with d in Eq. fl3~8]). It is 
important to stress that the contributions from new physics are generally different for the 
B s and E>d systems. The ratio C s /C d can therefore be sizeably different from unity and has 



6 



to be taken into account in model independent analyses of the unitary triangle when using 
Eq. (0). 

In the numerical analysis we will require X s d to lay in the interval (30 4- 60) in order to 
satisfy the lower bound AMb b > 14.9 ps _1 and have an observable B s -B s mass difference at 
the same time. 



3 An example: SUSY gluino contributions 



In order to present a definite new physics model which contributes sizeably to the CP 
asymmetries in B — > (j)K, we turn to a variant of the MSSM with complex off-diagonal squark 
mass terms. In particular we will focus on the following entries: (Tn LL ) 23 , { m LR,)23> { m RR,)23 



and 



m 



RL 1 23' 



The other parameters of the model are the common mass of the squarks m, 



and the gluino mass nig. Following a common practice, we consider the normalized insertions 



°23 



2\d 
123 



to fh 2 . 



given by the ratios of the various 

It is well known, and in the next section we will give a detailed quantitative analysis of 
the question, that E(B — > X s j) puts severe bounds of order O(10~ 2 ) on the LR and RL 
insertions while the impact on the LL and RR ones is rather mild. This strong result follows 
from the nig /nib chiral enhancement of the (823) lr,rl contributions to the Wilson coefficient 
CV 7i 8 9 . Note that this chiral factor is absent in the in the WC's that govern the B s -B s mass 
difference and the B — > <pK amplitude. For the latter observables, in fact, all the mass 
insertions enter with similar weight: therefore, the only insertions that can play a role are 
the LL and RR ones. Since the analyses for both insertions give the same results, we will 
consider explicitely only ((^ll- 

The expressions of the SUSY contributions to the various coefficients are [EHI] 



C 3 

c A 

C 8 g 
qVLL 



'23>LL 



2V2G F m 2 \V tb V t * s 



tb v ts\ 
23 J LL 



15 1 

gB^Xgg) ~ -B 2 (Xgg) - -P(X~gg) 
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2V2G F m 2 \V tb V t * s \ 

a 2 s (S d 23)LL 

2V2G F m 2 \V tb V t * s 



--B^x-gz) + ~B 2 (x M ) + Pi 



Xgq) 



a 2 (5 d 



tb v ts\ 
'23 )LL 



2V2G F rh 2 \V tb V ts] L 

nUsiS^LL 16 , 

V2G F rh 2 \V tb V t l\ 9 92{X ^ 



— B^Xgq) + —B 2 (Xgg) ~ -P{Xgg) 

2 7 

-~B X {X M ) + — B 2 (Xgg) + P(Xgg) 



f(x ) + ™f<x )) 

2|y v*\ 2 \9 9q ^ m ' ^2 9 J ' 



(39) 
(40) 
(41) 
(42) 
(43) 
(44) 
(45) 



7 



where x g g = m~/m? (we assume a common squark mass). The loop functions can be found 



in Refs. [21], 2TJ. Note that we have included the phase of V tb V t * (which is basically a minus 
sign) in the definition of the mass insertion. 



4 Numerical analysis 

We begin with the constraints imposed by the B — > X s 7 branching ratio and CP asymmetry. 
We assume the absence of so-called accidental cancellations between different large contribu- 
tions; therefore, we present the analysis assuming the presence of only one insertion per time 
besides the SM. If all of the insertions are substantial at the same time, the allowed ranges of 



some combination of them can be sizeably enlarged (see for instance Ref. |T2] where the case 
of multiple real mass insertions is considered). In Fig. [I] we take x^ = 1, rh = 250, 500 GeV 
and require \5\ < 1; we also show the impact of using the 68% C.L. and 95% C.L. constraints 
on the B — > Xg'j branching ratio. The corresponding bounds induced by the CP asymmetry 
are much weaker and do not impact, at the moment, the allowed areas: future experimental 
improvements can substantially modify this picture. We plot the resulting regions in the 
[Re5, Im5] plane. Note that the contours scale with rh 2 . As mentioned before, the bounds 
on the LL insertion are not very strong. Even for very light squarks (i.e. rh = 250 GeV) it 
is possible to largely evade the constraint if the imaginary part of the insertion is 0(1). This 
reflects the absence of interference with the (real) SM contribution. On the other hand, the 
corresponding bounds on the LR insertion are, as expected, of order O(10 -2 ); moreover, as 
it follows from the figure, if \Slr\ is larger than O(10 -3 ), a substantial correlation between 
real and imaginary parts is required. For the 5r R and 5 RL insertions the results are similar. 

Next we analyze the possible deviations of the CP asymmetries in B — > <pK from the 
SM expectation and the relation to the B s -B s mass difference which is also altered by the 
insertions. In Fig. |2| we plot the correlation between the B s -B s mass difference and SA mix = 
j^mix _ j^mix.sM _ ^y e scan over ^ ne j n p U t SUSY parameters in the ranges 

rh E [250, 1000] GeV (46) 

K^WIeM (47) 



arg(5 2 d 3 ) LL G[0,27r], (48) 



for different values of x^, the gluino-squark mass ratio. We require each point to satisfy the 
B — > Xg'-f 95% C.L. constraint and to give a B s -B s mass difference in the range 30 < X S( i < 
60. The various regions correspond to the limiting case: for a given ratio x gq -, no point lays 
outside them. For x gq - > 1 we find that the deviations from the SM expectation remain 
below 0.05. For smaller values of Xgq much larger and thus observable contributions are 
possible. This strong dependence on the ratio between the gluino and squark masses is due 
to the particular dependence of the loop functions on Xgq. Moreover, the presence of definite 
bands in the [X S d, 5A mix ] plane is due to an interplay between the B — > X^ constraint and 
the requirement of fixed AMb s - 

As previously stated, 5A mix measures the deviation of the mixing-induced CP asymmetry 
from the standard model prediction. Since all new physics in the Bd-Bd mixing amplitude 
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will affect the CP asymmetries in the decays B — > J/^/K and S — > in the same way, 
<5A mia; is equal to sin 2(/3 + 9 d + 9 a) — sin 2(/3 + 9d) where 9d can be generated by other SUSY 
couplings that we do not consider here and 9a is the phase of the B — > (pK decay amplitude. 
In other words, SA mix is the difference between the mixing-induced CP asymmetries in the 
decays B — ► J/^/K and B — > cf>K (see also the discussion in Sec. 0). 

The new weak phases in the B — > amplitude also leads to a non vanishing direct 
CP asymmetry, A d * r , which can be measured for instance in decays of charged B- mesons. 
This asymmetry depends crucially on the presence of a strong rescattering phase, provided 
by the term P of Eq. fllTD- Nevertheless, we find that the new physics contributions to 
the two asymmetries are strongly correlated. This can be understood as follows. Let us 
parametrize the B — > (pK decay amplitude as A\^ x e l5x + ^e*^ 2 e tS2 where fa and Si are 
weak and strong phases respectively; note that this parametrization is arbitrary but that 
all physical observables do not depend on its choice. In the SUSY model that we consider, 
the weak phases fa are entirely due to the imaginary part of the mass insertions. Using the 
above parametrization, we find that the ratio A dtr /SA mix has an extremely tiny dependence 
on the phases fa; therefore, the correspondence between A dir and SA m%x is almost one-to-one. 
In Fig. [3], we explicitely show this correlation. 
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Figure 1: Constraints on the complex mass insertions {S^ll an d (S^lr coming from 
BiB — > X s j). We consider ra-frh 2 = 1. The contours scale as fh 2 . In the plot we explicitly 
show the cases fh = 250, 500 GeV. 
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Figure 2: Correlation between the B s -B s mass difference and the possible deviations 5A mix of 
the mixing-induced CP asymmetry from the SM expectation. The various regions correspond 
to different values of x g g. 
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Figure 3: Correlation between 5A mtx and the direct CP asymmetry A 



5 Conclusions 



We considered the effects of sizeable flavour changing entries in the squark matrices on b — ► s 
transitions. In particular we allowed for complex values of the relevant off-diagonal elements. 

We first investigated the bounds that these entries must obey in order to satisfy the b — > 
57 data. Fig. [I] shows that the inclusion of complex values for the mass insertion parameters 
strongly enlarges their allowed regions. In fact, the absolute values of the insertions can 
be much larger than in the existing literature where only real couplings were considered. 
Moreover, there are interesting correlations between real and imaginary parts which may 
give important hints on the structure of the underlying theory. 

We then considered the influence of the new terms on the CP-violating asymmetry in the 
decays B — > <pK and on the B s -B s mixing (the Bd-Bd mixing is not affected by the terms 
we are interested in). In Fig. |2|, we plot the deviation SA mix of the mixing-induced CP 
asymmetry in B — > <pK from the SM expectation versus the ratio X s d- We see that 8A mix 
can reach the 20% level in some corners of the parameter space; on the other hand, deviations 
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of order 0(10%) are easily possible with moderately light squark and gluino masses. Note 
that such large contributions are possible only for configurations in which m g < nig. In a 
similar fashion, Fig. |3] shows that the direct CP asymmetry A dir can receive contributions 
of the same order of magnitude. We stress that in this framework effects on B — > J/^/K 
decays are expected to be tiny. 

The B s -B s AMb s is very sensitive to the mass insertions we consider, while AM# d remains 
unaffected. This implies that the determination of V t d from ratio AMbJ AM^ may be 
misleading. Moreover, as it follows from Fig. |2| and |3], an experimental determination of 
AMb s in excess with respect to the SM prediction together with sizeable SA mix and A dir 
would be strong signatures in favour of this kind of models. During the next year, the B- 
factories BABAR and BELLE will gather enough luminosity to study the CP asymmetries 
in B —>■ <pK decays and will test this class of SUSY models soon. 
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